Nonconservative forcing and diffusion in refractive optical traps 
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Refractive optical trapping forces can be nonconservative in the vicinity of a stable equilibrium 
point even in the absence of radiation pressure. We discuss how nonconservative 3D force fields, in 
the vicinity of an equilibrium point, reduce to circular forcing in a plane; a simple model of such 
forcing is the refractive trapping of a sphere by a four rays. We discuss in general the diffusion of an 
anisotropically trapped, circularly forced particle and obtain its spectrum of motion. Equipartition 
of potential energy holds even though the nonconservative flow does not follow equipotentials of the 
trap. We flnd that the dissipated nonconservative power is proportional to temperature, providing 
a mechanism for a runaway heating instability in traps. 

PACS numbers: 05.60.Cd, 37.10.Vz, 42.50.Wk, 45.10.Na, 05.40.Jc, 42.15.-i, 87.80.Cc 



INTRODUCTION 

Nonconservative fields of optical force on optically 
trapped particles have long been predicted to occur [T]. 
These force fields are locally nonpotential in character so 
that net work is done on a particle even in microscopic 
closed paths. This differs from, say, electromotive force 
around a circuit where only the global potential is not 
definable[2]. In a locally nonconservative force field, ex- 
ternal power is continuously coupled to particle motion, 
leading to dissipation and heating even when the particle 
is localized about a stable point of the force field. 

Radiation pressure on trapped particles can be one 
source of nonconservative forcing. In recent experiments 
[31 [3] nonconservative toroidal circulation of optically 
trapped particles caused by axial radiation pressure was 
observed and theoretically analyzed. Recent theoretical 
work j5] has shown how, for nonspherical objects, non- 
conservative motion arises in coordinates of angle and 
translation. 

Here we show that a simple spatial 3D nonconservative 
force field, circular rather than toroidal, occurs even in 
the absence of radiation pressure. As a physical model 
leading to this force, we describe a sphere trapped by 
refracting rays. Within any such circular-forcing model, 
and generalizing to an anisotropic trapping force, we de- 
rive the spectrum of motion and the thermal signatures 
of nonconservative circular forcing. 

GEOMETRY OF NONCONSERVATIVE FORCES 

Before narrowing the discussion to optical trapping, we 
ask what simple generic statements can be made about 
nonconservative forces. Powerful classification schemes 
are available from differential geometry [SHE] but to use 
these one must correctly identify forces as fields, not of 
vectors, but of differential 1-forms. This means that force 
components such as f^, fy, and fz properly take their 



meaning from the work differential, or work 1-form, 

UJ = f^dx + fydy + fzdz (1) 

which is the integrand for evaluating work along any cho- 
sen path. The 1-form uj is said to be exact if Eq. ([T]) 
equals the differential of some potential function, uj — 
— d^; otherwise uj is inexact. For our purposes an exact 
1-form is the same thing as a conservative force. 

In Euclidean space one can reinterpret f^, fy, and fz 
as components of a vector. But as a 1-form field, Eq. ([T]) 
is constrained by Darboux's theorem [6l [7] , which states 
that every 1-form field uj can be reduced by some choice of 
general coordinates {qi, q2, q^, ■ ■ .) to a, shortest canonical 
form UJ = uj^^"^ belonging to the sequence 

^('^ = dgi, (2) 
^ q^Aq2, (3) 
w(^) dqi + gadgs (4) 

In more than three dimensions the list continues with 
o;*^^^ = qidq2 + qzdq^, and so on. 

The case uj = uj'^^\ i.e. uj — dqi, is exact, i.e. conserva- 
tive. Not only is qi a coordinate, but —qi is the potential 
function for uj. The case uj = uj^'^\ or w = qidq2, is re- 
ducible to the exact case by an integrating factor (1/gi), 
since [l/qi)uj — dq2. In two dimensions this exhausts our 
list, showing that all 1-forms (i.e., differentials) in 2D are 
either exact or integrable, a well-known and useful fact 
in thermodynamics. 

In 3D, exact (conservative) and integrable force fields 
of types UJ = uj'^^\ and uj = w^^-' can still occur, but the 
most general possibility is the nonintegrable 1-form uj^'^\ 
i.e. UJ — dqi + q2dq^. The form of a generic force field is 
thus more restricted than Eq. ([T]) would suggest. In the 
case of optical trapping, the optical force field includes 
a part derivable from a potential, which we use for our 
first coordinate: qi = — $. Following Eq Q, we write 

uj = fd0- d$ (5) 
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FIG. 1: A nonconservative 3D force field near a stable equilib- 
rium point. The circulation of force may be viewed as purely 
circular in some plane, (a) Vector flow-field point of view, (b) 
Combed-hair-on-a-sphere point of view. 

SO that 53 = 6* is a coordinate which we will below take 
to be an angle, and q2 = /, the angular force conjugate 
to 6. The Frobenius theorem [51-0 states that in terms 
of the exterior product A and exterior derivative d, the 
work 1-form will be of the nonintegrable form of Eq. ([5| 
whenever uj A dcu = — d$ A d/ A d0 7^ 0. 

In the neighborhood of an equilibrium point, where all 
components of a force field vanish, it is natural to in- 
terpret 9 in Eq. ^ as an angular variable, motivated 
by the reduction — yda: -I- xdy = p^d9. Nonconservative 
forces near equilibrium are circulations aligned with some 
plane that contains the equilibrium point. The noncon- 
servative portion of the force can be taken as strictly 
circular near the equilibrium point in any coordinates, 
since a noncircular 2D force pattern in x and y can be 
written as 

— aydx + bxdy— + b)p^ d9 — ^{a — b) d{xy) (6) 

and the last term can be absorbed into — c?"I>. Thus the 
work 1-form ^ is actually of the type ^ , with 9 a cir- 
cular angle in terms of x and y. 

Thus, to study a nonconservative 3D force about an 
equilibrium point, it is sufficient to consider a circular 
pattern of force in some 2D plane containing the equilib- 
rium point. 

From another point of view, this picture is consistent 
with equilibrium-point analysis of a vector field f , viewed 



as flow towards a fixed point [S] (Fig. 1). Briefly, with 
coordinates qi we have fi « MijSqj near the fixed point, 
defining some real matrix M. The eigenvalues of M are 
roots of a cubic equation with real coefficients, and gener- 
ically yield one real and two complex conjugate roots 
whose real parts, together with their eigendirections, cor- 
respond here to a conservative force field. The imaginary 
part of the complex pair of eigenvalues defines circulation 
in some plane. The eigendirections of M may be stretched 
and skewed, corresponding to the non-orthogonal coor- 
dinates in Q. 

From yet another point of view, any vector field re- 
stricted to the surface of a sphere (about the equilibrium 
point) will circulate about the sphere in a 2D fashion, 
according to the familiar "combed hair on a sphere" the- 
orem of mathematics [8] (Fig. 1). 

NONCONSERVATIVE TRAPPED-SPHERE 
MODEL 

For a "toy model" realization of a circular nonconser- 
vative force we consider an optically trapped sphere much 
larger than than the wavelength of light, so that geomet- 
rical optics applies. In our simplified model we assume 
no reflections (perhaps due to a graduated-index bound- 
ary) and we trap the sphere with only a few rays passing 
almost centrally through the sphere, so that refraction 
and trapping force can be calculated in a paraxial-ray 
approximation. Our model is certainly contrived, but 
is very simple to analyze. In more realistic situations, 
whenever there is chirality (handedness) of a system of 
light rays near an equilibrium point, we still expect some 
degree of circular nonconservative forcing to occur. 

We need only consider rays within the plane of 
Fig. 2(a), containing the sphere center. In the figure the 
displacement of the sphere center from the ray is exag- 
gerated. The force of the refracted ray on the sphere due 
to momentum transfer is 

f = --(k'-k), (7) 
c 

where Iq is the power of the ray and k and k' are unit 
vectors for the original and refracted rays. If A9 is the 
angle between k and k', then 

f = -— [(cosA6i-l)k+(sinA6i)f_L] (8) 
c 

where fj^ is the direction of rj^ (see Fig. 2(a)). We set 
sin A9 « A9 and cos — 1 « — A0^/2, and use paraxial 
ray transfer matrices [10] to express A9 and f in terms 
of the displacement r of the sphere. These matrices act 
on a vector {A9, h) whose components are the angle with 
and displacement from the optical axis. If matrices Fjn 
and Fout describe refraction at the convex surfaces of our 
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FIG. 2: (a) Refraction of a single ray by a sphere. In the 
paraxial limit, rx is much smaller than the radius of the 
sphere. The force f acts to return the sphere center to the line 
of k, but also has a component along k. (b) System of four 
rays for nonconservative trapping. The limit b <^ a allows a 
paraxial approximation for each ray. 



sphere of radius a, and and matrix describes parax- 
ial propagation by a distance 2a, the ray matrix for the 
sphere is 



FoutP2aFin — 



1 - n/2 -(n - l)/a 
a l-n/2 



(9) 



with n the relative refractive index of the sphere. In our 
case, h = ~r± gives 



(10) 



for the angular deviation. For simplicity we set n = 2 to 
find 



2a2 



(11) 



Eq. ( 11 1 resolves f into a conservative restoring force, pro- 



portional to — rj^, and a nonconservative quadratic force, 
pushing the sphere along the ray direction k whenever 
the sphere center is displaced off the ray axis. 

The system of four rays shown in Fig. 2(b) results in a 
purely circular nonconservative force on the sphere. De- 
spite the clockwise circulation of rays, there is no torque 



on the sphere within this ray-optical model: it is the 
nonconservative net force on the sphere that we study 
here. We separate oppositely directed ray pairs by a dis- 
tance 25, exaggerated in the figure, with 6 ^ a so that 
our paraxial analysis remains valid. The ray in the -l-a: 
direction is described by 

k = x, r± = {y - b)y + zz (12) 

and similarly for the other rays. Summing the forces from 



Eq. (11), the total force on the sphere is 
'b 



f = 



2/n 



-{xy - yx) - (xx + yy + 2zz) 



The work differential, or work 1-form, is 



— - p dt* - d 



(13) 



(14) 



which is of the form of Eq. ([s]), w = fdO — d"I>. 



NONCONSERVATIVE MOTION OF A 
DIFFUSING PARTICLE 

To look for spectral features of trapping in a noncon- 
servative field, we consider diffusion through a fluid of a 
particle with drag coefficient 7. At temperature T, the 
diffusion coefficient is ^/k^T^ an Einstein relation |11) . 
Near an equilibrium point we consider 2D motion r(i) in 
a plane of circulation where the nonconservative compo- 
nent of force is 



f = ^zxr = 



(15) 



where ^, the strength of the circular forcing, has units of 
force per distance. As a 1-form, 



(16) 



First we check that there is no pathology regarding work 
done by the nonconservative force at small length scales, 
where a diffusing particle executes spatial cycles with a 
diverging frequency. For this, imagine the particle con- 
strained to move on a circle of radius e, like a bead sliding 
on a circular wire. In the absence of the external force, 
the time to diffuse around the circle is 5t ~ {'j/k^T) e^. 
In the presence of a nonconservative force, the rates 
of positive and negative circling are enhanced and sup- 
pressed according to 



k^T 

— 2 



± 



(17) 



where AMK is the work done by the nonconservative force 
in a cycle; for our force this is /S.W = 2nef ~ ^e^. 
The net rate is /S.v = v^ — ^ ISW j 76^ , implying 
that the dissipated power is 



(18) 
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We conclude that the power supplied by the nonconser- 
vative force vanishes at small scales, as the square of the 
spatial scale of the motion. 

We turn to the spectral properties of the diffusing par- 
ticle, and first establish formalism pjj with force and 
motion only in the x direction. For a force f{t) and mo- 
tion x{t) that extend over all time, the spectral density 
{fx)^ is defined by |l2j 



{f^xl,)=2'K{Jx)^5{u:~u:') 



(19) 



where fi^ and Fourier transforms of functions trun- 

cated outside of a time T. Throughout, (. . . )(^ specifies 
a spectral density, whereas (...) is an average. Standard 
manipulations |12j . or more simply the formal replace- 
ment ifujx'^,) = T {fx)uj, give the rate of work done by 
/ as 



/dW\ _ 1 
\^/ ~ 27r 



1 

2n 



iuifx)^ dio. (20) 



Spectral densities integrate to correlations or mean 
squared fluctuations [H] In the ID case, if a thermal 
Nyquist force N{t) acts on the object, with {N^)ui = 
2^kj^T (with 7 the drag coefficient) then in an external 
conservative force / = —k^x, with — Hx/l, the equa- 
tion of motion N — Kj.x + "/x gives us 

= ~ij{uj + iax)x^ (21) 
{N^N:) = j\al+uj^){x^x:} (22) 

The replacements (N^N*,) = T {N"^)^ and {x^x*^,) = 
T (x'^)uj immediately lead to 



2, _ 2KTh 



This Lorentzian spectrum integrates to 



(23) 



(24) 



as it should by the equipartition theorem [131. Since 
{fx)u — ~Kx{x'^)uj is an even function of lu, the rate of 
external work vanishes. 

Now include a circular external force. With an 
isotropic conservative force — Kr, the vector equation of 
motion N = Kr — ^zxr + 7r gives us 



= {a - iiu - ri zx) 



(25) 



with "x" a cross product, and where a — k/7, and 
rj = ^/7. In the absence of thermal forces, the parti- 
cle spirals to the origin with angular frequency 77 and 
a radius proportional to exp(— at), similar to Fig. 1(a). 
Using the column vectors 











1 ^UJ 





(26) 



we can analyze an anisotropic confining potential, with 



Eq. (21) generalized to 



—irj Lj+iay 



(27) 



and ay = Hy/j- To find the power spectrum we form the 
averaged product 



rt \Mt 



(28) 



a generalization of Eq. ( 22 1 . In terms of spectral densities 
we have 



27fcBTl = 72M(rr^)„M^ 



(29) 



where 1 is the identity matrix, and the spectral density 
matrix is 



(rrT) 



(x^)^ {xy)^ 



(30) 



Eq. (29) is easily solved for (rr^)^^, giving diagonal ele- 
ments 



7 [uj^ + {al+a'^ — 2ri'^)uj'^ + {ri'^+axay)'^] 



(31) 

which generalizes Eq. ( 23 ); {y'^)uj is the same with ax and 
ay interchanged. 



Eq. (31 ) represents the power spectrum of a;-motion of 
a trapped particle with a nonconservative force (77 7^ 0) 
and 2D anisotropy (a^ ay) both in the xy plane. 
Nonconservative forces make this power spectrum non- 
Lorentzian even in an isotropic trap. Other interest- 
ing effects such as inertial hydrodynamics and material 
properties [Mj [15] will of course also invalidate a sim- 
ple Lorentzian spectrum. Carrying out the uj integral in 



Eq. ( 24 ) , we obtain 



2 , ^2 

Ky 



(32) 



For (y^), Kx and Ky are interchanged. In the isotropic 
the nonconservative circulation pre- 
serves {x"^) — k^T I n. For large values of ^, 



2k^T 



(33) 



and the potential is effectively being averaged by rapid 
circulation. However, for all values of ^, Kx, and Ky we 
find 



(34) 



so that the equipartition of potential energy is unaf- 
fected by nonconservative circulation, even though the 
nonconservative flows do not follow equipotentials in an 
anisotropic trap. 
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Details of the flow pattern follow from the off-diagonal 
element of (rr^)^^, 



r]{ay-aa:)-2ir]uj 



7 [uj'^ + {al+a'^ — 2ri'^)uj^ + {ri^+axayy 



(35) 

obtained by solving Eq. (29). Integrating this over w 
yields 



(36) 



Using this result we can generalize Eq. ( 32 1 to mean 
squared radius at any 2D angle. For Gaussian statistics 
such as we have here, the probability distribution 



P{x,y) cx exp[— ir^Gr] 



(37) 



yields a matrix of averages (rr^) = G \ and a unit ellipse 
defined by 



1. 



(38) 



As a function of angular direction 0, the squared radius 
of the ellipse is 



{x'){y')-{xyf 

(y-^) cos^ 9 — 2 (xy) cos 9 sin 9 - 



ix"^) sin^ ( 



into which (a;^), (y^), and {xy) can be inserted from 
Eqs. (32) and (|36|. The resulting ellipses are plotted 



in Fig. 3 for various values of ^. The r.m.s. radius drawn 
in the figure is necessarily a line of flow. For increasing 
5 the flow pattern does not follow the ^ = equipoten- 
tial lines, but tilts by 0^ in the direction of the circular 
forcing, with 



tan26'f 



2e 



(40) 



The flow pattern also decreases in eccentricity with ^ 
until, at high circular forcing, we recover the circularly 



symmetric distribution described by Eq. 33 



The rate at which work is done on the system is 



(41) 



where (f • v)^ is a spectral density. In terms of Fourier 
transforms, 



Ky/ — 5 










1 = 00 



FIG. 3: R.m.s. radius (r^)^^^ (solid curves) of the position 
distribution P{x, y) of a diffusing particle in a trap with stiff- 
ness anisotropy Ky/Kx = 5, for various values ^ of counter- 
clockwise circular forcing in the xy plane. The solid curves are 
lines of fiow which, for ^ 7^ 0, do not follow the ^ = equipo- 
tential contour (dotted curves). Circular flow at ^ = 00 (final 



(39) panel, solid curve) has the r.m.s. radius [IksT / {Kx+K,y\ 



11/2 



Using Eq. ( 35 ) for the right-hand side and performing 



the integral in Eq. (41) gives 



dW\ _ 4A:bT^2 



(44) 



We find that the rate of work dissipated in the system by 
the circular nonconservative force is directly proportional 
to temperature, reflecting the enhancement of induced 
drift velocities by thermal spreading from the zero-force 
point. In contrast, nonconservative toroidal circulation 
due to axial radiation pressure O H] leads to a squared 
dependence on temperature. Depending on the rate of 
outward heat flow from the optical focus [IS] the pro- 
portionality of dissipation to temperature raises the pos- 
sibility of an instability towards runaway heating by a 
nonconservative component of the trapping force. 



SUMMARY 





T 


■ 1' 






yuj 




-1 







= ioj^ix^yl - yu:X*J 
= -2^uj\mx^y^ 

from which we infer the spectral density 

(f • v)„ = ~2S,uj Im (xy)^. 



(42) 



(43) 



We have characterized the simplest 3D nonconserva- 
tive force field, near a stable fixed point, as an anisotropic 
conservative force plus a nonconservative circular force. 
We constructed a simple trapping model exhibiting such 
forcing, and presented signatures of particle displacement 
and nonconservative dissipation in the presence of non- 
conservative forcing. In most optical trapping experi- 
ments, simple circular forcing about a stationary point 
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may not be easily observable given the more dominant ra- 
diation pressure effects that have been observed by others 
[31 11] . Specifically designing chiral trapping geometries 
for particles and minimizing reflections may allow obser- 
vation of circular forcing. More generally, however, the 
geometrical picture we have constructed and the ther- 
mal and spectral results we have obtained may be useful 
in other situations, even outside of optics, where locally 
nonconservative microscopic forces act. 
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